High prevalence of pre-treatment and acquired HIV-1 drug resistance mutations among non-citizens living with HIV in Botswana.

Background: Approximately 30,000 non-citizens are living with HIV in Botswana, all of whom as of 2020 are eligible to receive free antiretroviral treatment (ART) within the country. We assessed the prevalence of HIV-1 mutational profiles [pre-treatment drug resistance (PDR) and acquired drug resistance (ADR)] among treatment-experienced (TE) and treatment-naïve (TN) non-citizens living with HIV in Botswana. Methods: A total of 152 non-citizens living with HIV were enrolled from a migrant HIV clinic at Independence Surgery, a private practice in Botswana from 2019-2021. Viral RNA isolated from plasma samples were genotyped for HIV drug resistance (HIVDR) using Sanger sequencing. Major known HIV drug resistance mutations (DRMs) in the pol region were determined using the Stanford HIV Drug Resistance Database. The proportions of HIV DRMs amongst TE and TN non-citizens were estimated with 95% confidence intervals (95% CI) and compared between the two groups. Results: A total of 60/152 (39.5%) participants had a detectable viral load (VL) >40 copies/mL and these were included in the subsequent analyses. The median age at enrollment was 43 years (Q1, Q3: 38-48). Among individuals with VL > 40 copies/mL, 60% (36/60) were treatment-experienced with 53% (19/36) of them on Atripla. Genotyping had a 62% (37/60) success rate - 24 were TE, and 13 were TN. A total of 29 participants (78.4, 95% CI: 0.12-0.35) had major HIV DRMs, including at least one non-nucleoside reverse transcriptase inhibitor (NNRTI) associated DRM. In TE individuals, ADR to any antiretroviral drug was 83.3% (20/24), while for PDR was 69.2% (9/13). The most frequent DRMs were nucleoside reverse transcriptase inhibitors (NRTIs) M184V (62.1%, 18/29), NNRTIs V106M (41.4%, 12/29), and K103N (34.4%, 10/29). No integrase strand transfer inhibitor-associated DRMs were reported. Conclusion: We report high rates of PDR and ADR in ART-experienced and ART-naïve non-citizens, respectively, in Botswana. Given the uncertainty of time of HIV acquisition and treatment adherence levels in this population, routine HIV-1C VL monitoring coupled with HIVDR genotyping is crucial for long-term ART success.

Predicted resistance to broadly neutralizing antibodies (bnAbs) and associated HIV-1 envelope characteristics among seroconverting adults in Botswana.

We used HIV-1C sequences to predict (in silico) resistance to 33 known broadly neutralizing antibodies (bnAbs) and evaluate the different HIV-1 Env characteristics that may affect virus neutralization. We analyzed proviral sequences from adults with documented HIV-1 seroconversion (N = 140) in Botswana (2013-2018). HIV-1 env sequences were used to predict bnAb resistance using bNAb-ReP, to determine the number of potential N-linked glycosylation sites (PNGS) and evaluate Env variable region characteristics (VC). We also assessed the presence of signature mutations that may affect bnAb sensitivity in vitro. We observe varied results for predicted bnAb resistance among our cohort. 3BNC117 showed high predicted resistance (72%) compared to intermediate levels of resistance to VRC01 (57%). We predict low resistance to PGDM100 and 10-1074 and no resistance to 4E10. No difference was observed in the frequency of PNGS by bNAb susceptibility patterns except for higher number of PNGs in V3 bnAb resistant strains. Associations of VC were observed for V1, V4 and V5 loop length and net charge. We also observed few mutations that have been reported to confer bnAb resistance in vitro. Our results support use of sequence data and machine learning tools to predict the best bnAbs to use within populations.

Predicted broadly neutralizing antibody (bnAb) resistance and associated envelope characteristics of adults with HIV-1 seroconversion in Botswana.

We used HIV-1C sequences to predict (in silico) resistance to 33 known broadly neutralizing antibodies (bNAbs) and evaluate the different HIV-1 env characteristics that may affect virus neutralization. We analyzed proviral sequences from adults with documented HIV-1 seroconversion (N=140) in Botswana (2013-2018). HIV-1 env sequences were used to predict bnAb resistance using bNAb-ReP, to determine the number of potential N-linked glycosylation sites (PNGS) and evaluate env variable region characteristics (VC). We also assessed the presence of signature mutations that may affect bnAb sensitivity in vitro. We observe varied results for predicted bnAb resistance among our cohort. 3BNC117 showed high predicted resistance (72%) compared to intermediate levels of resistance to VRC01 (57%). We predict low resistance to PGDM100 and 10-1074 and no resistance to 4E10. No difference was observed in the frequency of PNGS by bNAb susceptibility patterns except for higher number of PNGs in V3 bnAb resistant strains. Associations of VC were observed for V1, V4 and V5 loop length and net charge. We also observed few mutations that have been reported to confer bnAb resistance in vitro. Our results support use of sequence data and machine learning tools to predict the best bnAbs to use within populations.

High prevalence of hepatitis delta virus among people with hepatitis B virus and HIV coinfection in Botswana.

BACKGROUND: Approximately 15-20 million people worldwide are infected with hepatitis delta virus (HDV), which is approximately 5 % of people with chronic hepatitis B virus (HBV). Sub-Saharan Africa has high HDV prevalence, leading to worse clinical outcomes among people who are HIV/HBV/HDV tri-infected. There are limited data on HDV prevalence among people with HIV (PWH) who are HBV-infected and uninfected in Botswana. We, therefore, determined HDV prevalence among PWH in Botswana. METHODS: This was a retrospective cross-sectional study utilizing archived plasma samples from PWH with results for HBV markers such as hepatitis B surface antigen (HBsAg), hepatitis B core antibody (anti-HBc), immunoglobulin M antibody to hepatitis B core antigen (IgM anti-HBc) and hepatitis B e antigen (HBeAg). Samples were categorized according to their HBsAg status and screened for anti-HDV antibodies. Total nucleic acid was extracted from samples with a single positive anti-HDV result, and HDV ribonucleic acid (RNA) load was quantified using the Altona Diagnostic RealStar® HDV RT-PCR kit. Statistical analysis was performed using STATA version 14.0 where p-values < 0.05 were considered statistically significant. RESULTS: The study cohort (n = 478) included both HBsAg positive (44 %) and negative (56 %) participants, with a median age of 42 [IQR; 41-43]. Anti-HDV prevalence of (15/211) [7.1 %, 95 % CI: 4.4 - 11.4] was recorded among HBsAg positive participants, all of whom were IgM anti-HBc negative, while 5/6 participants were HBeAg negative. HDV RNA load was detected in 11/12 (92 %) anti-HDV-positive participants. No HDV prevalence was recorded among participants who were HBsAg negative, therefore, the overall HDV prevalence was (15/478) [3.1 %, 95 % CI: 1.9 - 5.1]. HIV viral load suppression was statistically insignificant, irrespective of HDV status. CONCLUSIONS: We report high HDV prevalence among HBsAg-positive PWH in Botswana. Most HDV-positive participants had active HDV infection, therefore, we recommend HDV screening in this cohort to guide their clinical care.